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Abstract. The magnetic susceptibility of the one-dimensional Hubbard model 
with open boundary conditions at arbitrary filling is obtained from field theory 
at low temperatures and small magnetic fields, including leading and next- 
leading orders. Logarithmic contributions to the bulk part are identified as 
well as algebraic-logarithmic divergences in the boundary contribution. As a 
manifestation of spin-charge separation, the result for the boundary part at low 
energies turns out to be independent of filling and interaction strength and 
identical to the result for the Heisenberg model. For the bulk part at zero 
temperature, the scale in the logarithms is determined exactly from the Bethe 
ansatz. At finite temperature, the susceptibility profile as well as the Friedel 
oscillations in the magnetisation are obtained numerically from the density-matrix 
renormalisation group applied to transfer matrices. Agreement is found with an 
exact asymptotic expansion of the relevant correlation function. 

PACS numbers: 71.10.Fd,71.10.Pm,02.30.Ik 

1. Introduction 

The one-dimensional Hubbard model plays a central role in the understanding of 
interacting electrons in one dimension. The Hamiltonian 

L-l L ^ L 

arises naturally in the tight-binding approximation of electrons on a chain with L 
sites. In (1), the magnetic field h couples to the ^-component of the total spin. 
The interaction parameter u = U/\t\ > is the ratio between the on-site Coulomb 
repulsion U and the hopping amplitude t. Note that the eigenvalues of H are invariant 
under a sign change t ^ —t [1]. Furthermore, H is invariant under reversal of all 
spins and under a particle-hole transformation (the so-called Shiba transformation) 
[1]. Therefore we restrict ourselves here to positive magnetisation and lattice filling 
less than or equal to one. 
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The appealing simplicity of the Hamiltonian, combined with its Bethe-Ansatz 
solvability and its adequateness for field-theoretical studies are the reasons for its 
wide popularity. 

Recent experimental achievements in two areas additionally motivate our studies: 
On the one hand, the fabrication and characterisation of carbon nanotubes have shown 
that these can be considered as realisations of one-dimensional Hubbard models [2]. 
Especially, kinks in these quantum wires have been realised experimentally [3, 4]. 
A kink locally weakens the hopping amplitude at one specific lattice site in the 
Hamiltonian. Such a modification is known to be a relevant perturbation, which, 
at u > 0, is governed by a fixed point with zero conductance through the kink [5]. 
Thus at low energies the chain is effectively cut into two pieces. 

On the other hand, ultracold fermionic atoms in optical lattices can be described 
by a one-band Hubbard model [6]. Given the recent progress in realizing quasi one- 
dimensional bosonic quantum gases [7], it is likely that similar experimental progress 
will be made with fermions. 

In order to model these situations, we consider open boundary conditions in 
(1). Compared to the case with periodic boundary conditions, an additional surface 
contribution fs to the free energy occurs, defined as 

fB = hm (Fobc - Fphc) , (2) 

L — >QO 

where i^obc {Fpbc) is the total free energy for open (periodic) boundary conditions. In 
this work, we will focus on the magnetic susceptibility per lattice site x — Xbuik+Xs /L 
both atT = 0,/i>OandT>0,ft. = 0, at arbitrary fillings. 

For half filling, the bulk contribution Xbnik at T = ft = has first been given 
by Takahashi [8], where also the existence of logarithmic contributions at finite h is 
mentioned. Later on, Shiba [9] calculated Xbuik at T = /i = for general filling. 
The free energy at finite temperatures was given by Takahashi (for an overview and 
original literature, cf. the book [10]), and later by Kliimper (the book [11] contains a 
detailed account of this work) . However, it seems as if the explicit behaviour of Xbuik 
at T = 0, /i > and T" > 0, h = has not been derived so far. In this work, this gap 
will be filled. 

Although the boundary quantity Xb is only an C'(l)-correction to the total bulk 
contribution, it may become important in experiments if it shows a divergency with 
respect to the temperature T or magnetic field h. Indeed, such divergences have been 
discovered and analysed in the spin-1/2 Hciscnberg chain [12, 13, 14, 15, 16]. Since 
the isotropic spin-1/2 Heisenberg chain is obtained from (1) in the limit w ^ oo, 
related divergences are also expected to occur in (1). In the case of half-filling for 
T = 0, it has in fact been shown in [14, 17] that the boimdary magnetic susceptibility 
is divergent, Xb ~ l/{h\n^h), for ft — » 0. Exactly the same result for xb has also 
been obtained for the supersymmetric t — J model [18]. The OBCs do not only lead 
to 1/L-corrections but also brc;ak translational iiivariance. Therefore quantities like 
the magnetisation or the density become position dependent. Local measurements of 
such quantities then provide a way to obtain information about the impurity making 
theoretical predictions about the behaviour of such one-point correlation functions 
desirable. 

In section 2 we give the functional forms of both the bulk and the boundary 
contributions by using a field-theoretical argument. Leading and next- leading 
logarithmic contributions to the finite bulk susceptibility are found both at finite 
T and finite ft. On the other hand, the boundary contribution diverges in a Curie-like 
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way with logarithmic terms, where again we give both the leading and next-leading 
divergences. These results depend each on two constants, which are the spin velocity 
and the scale involved in the logarithms. At magnetic fields or temperatures much 
smaller than this scale the result for the boundary susceptibility becomes independent 
of the spin velocity and the scale and therefore completely universal. This is a 
manifestation of spin-charge separation as will become clear in the next section. For 
the bulk susceptibility, on the other hand, only the functional dependence on field 
or temperature will be universal for low energies. The value at zero field and zero 
temperature, however, is a non-universal constant which does depend on filling and 
interaction strength via the spin velocity. 

The spin velocity has been determined previously from the Bethc ansatz solution 
[9]. In section 3, for T = 0, the scale in the logarithms will be determined exactly as 
well. The calculation of boundary effects at T > based on the Bethe ansatz solution 
still remains an open problem, as for all Yang-Baxter intcgrable models (reasons for 
that are given in [15, 19] for the special case of the spin-1/2 XXZ chain). In section 
4 we therefore calculate the susceptibility profile xi^^T) and magnetisation profile 
s{x,T) in the asymptotic low-energy limit (that is. for large distances and small 
temperatures) by making use of conformal invariance. Due to the open boundaries, 
s{x,T) shows the characteristic Priedel oscillations [20, 21, 22]. To test the field- 
theory predictions, we perform numerical calculations in the framework of the density 
matrix renormalisation group applied to transfer matrices, which is particularly suited 
for impurity and open-boundary models. In the last section we will present our 
conclusions and discuss in which experimental situations the calculated boundary 
effects might become important. 

2. The low-energy effective Hamiltonian 

First, we briefly review the effective field theory for the Hubbard model following in 
large parts Refs. [23, 11]. From the effective Hamiltonian we then obtain the magnetic 
susceptibility at small magnetic field and low temperature. 

Let a be the lattice spacing. We introduce fermionic fields i){x) in the continuum 



where x = j • a and the usual splitting into left- and right-moving parts has 
been performed. The Fermi momentum depends on both the density n and the 
magnetisation s as A;ft = 7r(n-|- 2s)/(2o), kpi = 7r(n — 2s)/(2a) (the magnetisation is 
defined as s = (m| — mj^)/2, with = Ma/L being the density of particles with spin 
a). In the following, we will consider the zero-field case where kpi = kpi = 7rn/(2a). 
Eq. (3) allows it to introduce a Hamiltonian density W(a;), such that H = J Hdx. In 
terms of the right- and left-movers in (3) the kinetic part of the Hamiltonian (1) in a 
lowest order expansion in a becomes 



by 



y/d,ip{x) = ^/^{e''""'''R„{x) + e 



(3) 




(4) 



(7 



and the interaction part 



Hint = 4ua {: {r\r^ + (r\Ri + l\Li 



) : - : R\RiL\L^ : - : RIR^L\L^ : 




(5) 
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Here ":" denotes normal ordering. For brevity, the x-dependence of the operators has 
been dropped. The Fermi velocity is given by vp := 2asm(kpa). The second term in 
(5) represents backward scattering processes whereas the last term is due to Umklapp 
scattering. Only in the half- filled case, where kp = n/{2a), is the Umklapp term 
non-oscillating and has to be kept in the low-energy effective theory. For all other 
fillings it can be dropped. 

To make the Hamiltonian manifestly SU{2) invariant under spin-rotations one 
can also express H in terms of the following currents [24, 11, 25]: 



(7 

J° = - ^ : Rla^^Rf) : 



a,/3 



a, 13 



The free part (4) then reads 



Ho 



Vp 



10 



J' 



2-K 



: J J 



J J 



(6) 



As far as the interaction part (5) is concerned, we first consider the case n ^ 1, 
that is away from half-filling. Then Umklapp scattering can be ignored leading to 



Hi: 



J' 



■.T 



J J 



J J 



Ac : JJ : +A : J • J : 



(7) 



and coupling constants gc = ua, gg = — 4ua/3, Ac = 2ua and A = —Sua. 

Taking Eqs. (6), (7) together, we see that the Hamiltonian is a sum of two 
terms: one depending on the scalar currents J, J only (corresponding to charge 
excitations) and the second depending on the vector currents J, J (associated with 
spin excitations). 



He 



TTVp 



+ J 



: +Ac : J J 



(^+fl«):[j-J + J-j]:+A:J-J:. 



(8) 
(9) 



The charge and spin parts commute, [HcHs] = 0. 

Let us first focus on He- The term with coefficient gc gives rise to a 
renormalisation vp, yielding the charge velocity 

Vc = vp + 2ua/Tr. (10) 

Upon bosonising, the charge currents are written as [23] 



1 



(n + d.ip) , J 



1 



(n - d,^) 



V47r V47r 
where the boson field (p and the corresponding momentum H satisfy the canonical 
commutation relation [(/^(a;), n(a;')] = id {x — x'). Then 



n2 1 1 



^) + {d^pf (i + A 



By scaling (f' = <^\/^) n' = n/\/Kc, 

__ A^ 2ua 



(11) 
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this Hamiltonian is written as Wc = x + i'^x'p')'^ , which describes noninteracting 
fields. Note, that in this field-theoretical approach the Luttinger parameter Kc is 

calculated only up to the linear order in u. The same is true for Vc^s- Contributions 
in higher tt-order would occur if the perturbational approach is pursued further. 
Fortunately, the Bethe-ansatz solvability of the model allows it to calculate Vc,s,Kc 
exactly [11]. We will come back to this point in the next section. 

In Tis , the gs-tevm leads to a renormalisation of the spin velocity 

Vs = vf — 2ua/7r. (13) 

The remaining interaction of vector currents is a marginal perturbation. By setting up 
the corresponding renormalisation group equations, it turns out that it is marginally 
irrelevant (relevant) for sgn(A) < (sgn(A) > 0) [23, 26]. In our case, A = —Sua < 0. 
The remarkable point about this is that the spin part of the low-energy effective 
Hubbard model is identical to the corresponding expression for the XXX-Heisenbcrg 
chain [27, 28], whereas the charge part is described by free bosons (away from half 
filling). For this case, field theory has been employed [27, 28, 26] to obtain the bulk 
contribution to the magnetic susceptibility in the form 

/ 1 lnln£;o/£' IE \ 

» = ('^' 

where E = h,T (magnetic field or temperature), Eq = Hq^Tq is some scale and 
Xo := x{T = 0, /i = 0) is given by the inverse of the spin velocity. For the open 
XXX-chain, the boundary contributions have been found to be [12, 13, 16]: 

1 / InlnTo^^Vr 

''''^ ^ 12TlnT(^Vr V 21nT(^VT ^ InT^^V^ ^"') ^ ' 

1 / lnln/if,^V/i ll^'' \ , ^ 

XB{h)= , 1 + + . . . 17 

From the above considerations we conclude that the bulk and boundary contributions 
to the magnetic susceptibility in the Hubbard model are also of the form (14)-(17), 
where, compared to the XXX-model, ^o, ^O) lT,h, ItI renormalised by the 
charge part. 

Most interestingly, the pre- factor of xs remains unaffected by the charge channel. 
The divergent boundary contribution for T <^ or /i <C /iq is therefore completely 
universal 

1 /, InllnTl 



This can be understood as follows: As the charge- and spin-part of the Hamiltonian 
separate at low energies, the only way how the charge channel can effect the spin 
channel is by a renormalisation of Vs and Kg. The Luttinger parameter Kg in the 
spin sector, however, is fixed to Kg ^ 1 due to the SU (2) symmetry and cannot 
renormalise. The explicit calculations of the pre-factor of xb for the XXX-model in 
[12, 13, 16] show, on the other hand, that it does not depend on the spin- velocity Vg- 
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It therefore remains completely independent of the filling factor and the interaction 
strength. 

Let us now comment on the scales involved in Eqs. (14), (16), (17). Including 
the order 0{lii~^ E) in Eq. (14), this equation can be written as 

Xbuik(i?) _ ^0 ^1 _ - + + . . .J . (20) 

It is convenient to define a new scale Eq = Eq^'^'^ . Then, again up to the order 
O (In"^ E) , we have 

( 1 lnlni?o/£^ \ 
XhME) = XO 1 + + . . . , 

where the term ~ ln~^ E has been absorbed in the term ~ In^^ E by the redefinition 
of the scale. This procedure fixes the scale uniquely [26]. 

One proceeds analogously with Eqs. (16), (17) and obtains 

^"^"^^ - 12Tlnr(-)/r [' - 21nff)/T ^ " ' " j ^''^ 

XB{h) = 1 ^ °. ' +■■■] ■ 23 

Let us now turn to the half- filled case n = 1. The additional Umklapp term 
in (5) can also be bosonised and is proportional to cos \/8n(p. When we now again 
rcscalc the field (p' = (p\/Kc wc obtain a relevant interaction ^ cos y^Sn/Kcip' for 
any finite m > because Kc < 1 in this case. This means that the charge sector will 
be massive. Indeed, at u — > oo, the excitations of the Hubbard model at half filling 
arc exactly those of the XXX-chain [11]. The formulas (14,16,17) remain valid at 
half-filling as well. The leading term in Eq. (17), including the constant ho, was given 
in [17]. There, a phenomenological argument was found that generalises this result to 
arbitrary filling. The constant ho was left undetermined in the arbitrary filling case. 



3. Bethe ansatz 



In the framework of the Bethe ansatz (BA) solution, the energy eigenvalues of (1) 
arc given in terms of certain quantum numbers, the Bethe roots. These roots have 
to be calculated from a sot of coupled algebraic equations. In the thermodynamic 
limit, these algebraic equations can be transformed into linear integral equations for 
the densities of roots, with the energy being given by an integral over these densities. 
In this section, the Bethe ansatz solution is first used to verify the small-coupling 
expressions for Vc,Vs,Kc, cf. Eqs. (10,12,13). 

Afterwards, we obtain the magnetic susceptibility at T = 0. Therefore, we 
first analyse the integral equations in the small-field limit, thereby determining the 
constants in Eqs. (14) (for T = 0). The pre- factor xo has been calculated by Shiba 
[9] . Our essential new results are twofold: On the one hand, the leading /i-dependence 
of the bulk-susccptibility is calculated exactly at small fields, including the scale, for 
arbitrary fillings. On the other hand, the result for the boundary susceptibility (17) 
is confirmed within the exact solution. However, due to cumbersome calculations, the 
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constant 7)^ ' in Eq. (17) is left undetermined here. Finally, we consider some special 
cases and present numerical results. 

In order to introduce our notation, we shortly state the main results of the BA 
solution. For any further details, the reader is referred to [11], which also contains 
an extensive list of the original literature. The BA solution for the one-dimensional 
Hubbard model with open boundary conditions has been found by Schulz [29], based 
on the coordinate Bethe ansatz. The algebraic BA for this model has been performed 
in [30]. The BA equations read 



Ml 

■'■■^ Ai — sin fc.- 4- Ml. Xr 4- .sin fc.- — 11/. 



g2ife, (i+i) ^ 77 A; - sin kj -iu\i+ sin kj + iu 



Xi — sin kj + iu A; + sin ki — iu 

1=1 

^-r A; — siu /uj- — i(/ A/ + siu/ij- — iu -p-i- A; — A,„ — 21;/ A; + A,„ — 2iu 
■I--'- A; - sin A;,- + iu A; + sin kj + iu , ^, A; - A^ + 2iu A; + A^ + 2iu ' 

l = l,...,M^, (25) 
and the energy is given by 

JV 

E = -2^coskj. (26) 

Here, we analyse the groundstate, where the A'^-many fc^s and the A^f^-many A;s lie on 
one half of the real axis, except the origin. Although fc = 0, A = are solutions of the 
system (24), (25), they must not be counted in (26): These solutions correspond to 
zero-momentum excitations, which have to be excluded due to the broken translational 
invariance in the open system. | However, one can show that if kj,Xi solve (24), (25), 
then the same is true for —kj, —A;. One then "symmetrises" Eqs. (24), (25) by 
setting up equations for the sets {pi, . . . ,p2N+i} ■= {— fcjv, . • • , — fci, 0, fci, . . . , fcjv} 
and {vi,.. .,V2Mi+i} {-Am^, . . . , - Ai, 0, Ai, . . . , Am^}:5 



e 



2Mi+l 

2ip,(z,-n) Sinpj -Hm _ -pr vi - smpj - lu 



. , 3 = l,...,2N+l 27 

1)1 — sm r). 4- Mt 



svsipj — lu j"--^ vi — smpj -\- lu 



v, + 2m -p. .,-smp,+m^ ^ " ""^ = 1, 2M, + 1. (28) 

These equations can be solved analytically in the small coupling limit, cf. Appendix 
A. However, this solution has to be treated with care. It has been shown in [31] that 
a small-coupling expansion of the groimd-state energy in the thermodynamic limit has 
zero radius of convergence. This does not come as a surprise when considering again 
the low-energy effective theory presented in Sec. 2: At u = the interaction of vector 
currents in Eq. (9) changes from marginally relevant to marginally irrelevant. Thus 
in the following we perform the thermodynamic limit before considering any small- 
field or small-coupling approximations and compare with the results of Appendix A 
afterwards. 

X The wave function constructed from the coordinate BA [29] would vanish identically for fe = 0, A = 
0. 

■S fc = 0, A = axe included here to introduce homogeneous densities of roots. Their contribution is 
then subtracted later on, see Eqs. (32,33). 
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In the thermodynamic limit, one can set up equations equivalent to (27), (28), by 
introducing the density of ps, p{k), and the density of As, (t{v). These densities are 
solutions to the following set of coupled linear integral equations [29, 14] 

i-B 



11/1 \ 
p{k) = — \- T { cos A; ai (sin fc) +cosfc / ai{smk — v)a{v)dv (29) 

1 f'^ 

(^i"^) = YC'2{v) + / ai{v — smk)p{k)dk — / a2{v — w)a{'w)d'w, (30) 

L J-Q J-B 



where 

an{x) = nu/{'K{n'^u^ + x^)). (31) 

The integration boundaries are determined from the particle density n and the density 
of particles with spin down m^, 

" = ^ /.Q P(fc)dfc -±.,m^ = ^J^^ a{v)dv - ^. (32) 

Once these equations are solved, the energy density e is obtained from 

rQ 1 
■ / cosA;p(fc)dfc + -. (33) 
J-Q L 



e 

' -Q 



3.1. Velocities and Luttinger parameter 

Before proceeding further, we first make contact with the field-theoretical results 
(10,12,13) in the previous section. Since these concern only bulk-quantities, wc discard 
the 1/L-corrections in this subsection. We also set the lattice parameter a = 1 here. 

Furthermore, the results (10,12,13) have been obtained for densities n^l (such 
that both the charge and spin channels arc masslcss). Analogously, wc restrict 
ourselves here to densities away from half-filling. Within the BA, charge- and spin- 
velocities are calculated as 

1 i9Aec,s(A) 



decs dxec,sW 

Vc,s = - 



dpcs dxpc,s{X) 



\=B,Q 



(34) 

A=B,Q 



where ec.s is the energy of the lowest possible (i.e. at the Fermi surface) elementary 
charge/spin excitation and pc.s the corresponding momentum. After the second 
equality sign, these quantities arc parametrised by the spectral parameter X = k,v. 
Then pc,s is expressed by the densities p, a. The derivatives of the elementary 
excitations with respect to the spectral parameter are given by {e'^{k) = dkicik), 
e'siv) = dyCsiv)) 

.B 

e^(fc) = 2sinfc -|- cos fc / ai{siiik — v)e'g{v)dv (35) 

J-B 

rQ rB 
Ss(^) = / ai{v — smk)e'^{k)dk — / a2{v — w)e'g{w)dw. (36) 

J-Q J-B 

Furthermore, from Tic in Eq. (11) it follows that the Luttinger liquid parameter Kc is 
obtained from the charge susceptibility Xc at zero field, 

2K^ 
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Figure 1. Spin- (black point-dashed lines) and charge- (black full lines) velocities 
at fillings n = 1/4, 1/2, 3/4 (pairs from bottom to top). The dashed lines on the 
left are the low-ti results Eqs. (10), (13). The horizontal bars on the right indicate 
the limiting value Vc\u^<x = 2sin(7rn), [33]. Note that this asymptotic value is 
the same for both n = 1/4,3/4. By comparing Vc\u=o = = 2sin(7rn/2) with 
Vc\u^oo, one concludes that Vc(u) is maximal at a finite u for 2/3 < n < 1. 



The susceptibility Xc — d^n in turn can be expressed by e'(k)\k=B, p{B) and a 
related function [32]. Although the analytical solution of the integral equations in the 
limit M ^ is difficult to obtain due to singular integration kernels, these equations 
can be solved numerically to high accuracy. Figs. 1, 2 show the velocities and the 
Luttinger parameter at different fillings as a function of u, together with the analytical 
predictions (10,12,13) for small u. Agreement is found in all cases. 

Once Vs is calculated, xo is also known by virtue of Eq. (15). In the next section 
we describe how to obtain the /i-dependent corrections. 

3.2. Spin susceptibility 

An analytical solution of Eqs. (29), (30) is a challenging task due to the finite 
integration boundaries Q, B. To gain a first insight, consider the integral of Eq. (30) 
over the whole real axis, yielding cr(u)du = 2{n — mi) + 1/L and therefore the 
magnetisation 



Our aim is to perform a low-field expansion, i.e. an expansion around s = 0. From 
Eq. (38), this corresponds to an asymptotic expansion of a{v > B ^ 0). This 
expansion is done by generalising Shiba's approach [9], who calculated xo for pbc 
in a different way than via the spin velocity. 




(38) 
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Figure 2. Luttinger parameter Kc at fillings n = 1/4, 1/2, 3/4. The blue dashed 
lines are the low-u results Eq. (12). 



Substitute Eq. (29) into Eq. (30) to obtain 

'^W = ^(.ggiv) + SQ{v,0))+g^^\v) - |^ SQiv,v')aiv')dv' 

fQ 

Sq{v,v') := a2{v — v') — / ai{v — smk)ai{smk — v') coskdk. 

J-o 



Here is the v = 0-case of Qq' , defined by 



1 



(v) :— / — ai(w — sinfc) cos'' fcdfc. 

J-Q ^ 

Furthermore, 



SQ{v,v')a'^''\v')dv'. 



(39) 



(40) 



(41) 



Eq. (41) can be solved for <7^'^\ at the expense of introducing a new unknown function 
Mq{v,v'): 

poo 

(42) 



Mq{v, v) = Sq{v, v) - I Sq{v, v")Mq{v" , v')dv" 

J — oo 

'^^'^H") = 7(5?(«) + ^qKO)) + .g(;Hi;) 



Mq{v,v') 



j{9'-Q\v') + SQ{v'M + g'Q'W) 



dw'. 



(43) 



The one- dimensional Hubbard model with open ends 



11 



We consider now Mq{v, v')a{v')dv' , where (t{v') is given by the rhs of Eq. (39). 
Making use of Eq. (43) with v = 0, one obtains 



(j{v) = a^^'^iv) + [ MQ{v,v')a{v')dv'. 

J\v'\>B 



(44) 



Similarly, starting from ^^^a{v)(j^^'> {v)dv (where a'^'^^v) is given by Eq. (41)), one 



gets 



J\v\>B J-oo L L \ / 



-B 



(45) 



Let us now express the energy and the particle density in terms of these functions. By 
use of the definition (40), and recaUing Eqs. (32), (33), one arrives at 

I'B 

e- eo 



g^Q\v)a{v)dv + TT / g^Q\v)a^°\v)dv (46) 

-B J-oo 
" roo 

{v)a{v)dv - - / {v)<T^'^ {v)dv, (47) 



n — no = 



where eo, no are the energy and particle densities at i? = oo, i.e. zero magnetic field. 
Note that e, n are functions of both Q, B. First, we hold Q fixed so that both e and n 
change with varying B. At the end, we will account for the change in n and calculate 
the susceptibility at fixed n. 

Let us now distinguish between the bulk and boundary parts in the auxiliary 
function a^'^^ (an index b denotes the boundary contribution and should not be 
confused with the integration boundary B), 

1 



■ ^bulk ' j^"B 

Then with the help of Eq. (41) the second term in Eq. (46) is written as 



/oo /.c 
9^Q\v)a^'\v)dv= / 
-OO •/ — c 



JO), 



1 



1 



iv){l + -]+-SQ{v,0) 



(48) 



(49) 



The first term in Eq. (46) is reformulated with the aid of Eq. (45). Analogous 
calculations are done for Eq. (47). Putting everything together allows us to write 



e — eo = TT 



\v\>B 



a{v)a^'^\v)dv 



L 



\v\>B 



QqHv) + Sq{v, 0) CTbuik(u)dt; 



(50) 



n 



-no = ~ [ u{v)a^^\v)dv + ^ ( Wq\v) + SQ{v,Q)\a^,^,^{v)dv. (51) 

^ J\v\>B J\v\>B L J 



'|j;|>B J\v\>B 

We thus have to evaluate (j{v), cr^^'^^(v) asymptotically. Prom Eq. (41), 



oo /.Q 



ai(sinfc — v') cos A: 
Q 4u cosh 5^ (t" — sin k) 

^ cos" k 

.Q 4u cosh -^{v — sin k) 



a^''\v')dkdv' 



dk 



^ ai{smk) cosk 
Q 4u cosh ^ (w — sin k) 



dk 



(52) 
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has been defined. Eq. (52) can be solved for a^'^\ at the expense of introducing an 
unknown function LQ{t,t'), 

/sin Q 
K{t - t")LQ{t" ,t')dt" 
- sin Q 



sin Q 

sinQ y^sin Q 



dt / 

-sinQ J- 



dt'- 



SinQ J -SinQ 47ru cosh ^ (t; - i) 

1 + —) / cos" p K^^\sm.p — t')dp 
J-Q 

~T I 0,1 {sin p) COS pK^^\smp — t') dp + —K^^\t') 



(53) 



From the expression (53), the v ^ oo asymptotic behaviour of a^'^^ (v) can be read off, 
namely: 



(54) 



with an algebraic decay K^^^dti] oo) ~ l/(4w^). The quantity Iq^ is defined as 



-1 



Q 



cos^k 
2-K 



2L 



e5^r^'"*dfc- — / ai(fc) cosA;e5ii^'"*dfc 



Q 



/sinQ ^-sinQ i,/ / 1 \ f'^ 

dt ^e^'LQ{t,t') 1 + y / cos''p«W(sinp-t' 

-sinQ J -sinQ \ J-Q 



')dp 



1 TT 

— y / 7rai(sinp) cospK;(^^(sinp — f')dp + — 27rK(^^(t') 
L J-Q L 



(55) 



For later purposes, let us also separate this function into bulk and boundary parts. 



Q — -"Qibulk + ^ ^Q,B- 



From Eq. (44), we now calculate a(v) in the asymptotic limit. Therefore it is helpful 
first to reformulate Eq. (44) by writing 

oi(sinA: — w') cosfc 



Mq{v,v') = K^^\v-v') 

/•sinQ /"sin Q 



Q Au cosh ^{v — sin A:) 



dk 



dt dt'- -— -LQ{t,t') cos-'pK^^\sinp-t')dp. 

-SinQ i-sinQ StTM COSh ^ (u - f) J _q 

From this we conclude 

JUq{v + B,v' + B)+ J^q{v + B, -v' - B) « K^^\v - v') + k'^^\v + v' + 2B). (56) 
Therefore the equation 
g{v + B) = G^^\v + B) 

pOO 

+ / [Mq{v + B,v' + B)+ Mq{v + B, -v' - B)] a{v' + B)dv' (57) 
Jo 
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can be approximated by 

aiv + B) ie-ir("+B)j(o) + 1 + 5 



u 



+ 



f 

Jo 



1 1 



=: ^e-^(''+^)7^°)Pi(t; + B) + f PaC^^ + B) 

1j Zu 



-: e 



+ 



(u-?;') + «;^^'(^^ + ^'' + 2B) Pi{v')dv' 



P2{v) =: 2m(^)(v) 



{v - v') + k'^^\v + v' + 2B) P2{v')dv'. 



(58) 
(59) 
(60) 



Let us estimate the error involved in the above approximations. We will see later 
that B ^ — In ft, for small magnetic fields h. Corrections to the first term in Eq. (54) 
are higher-order exponentials and thus would add terms ~ /i^" to the susceptibility 
(the term taken into account here yields a constant contribution /i°). In the second 
term, higher-order algebraic terms have been dropped. These would contribute in 
order ~ l/(/iln"/i), n > 3, to the boundary susceptibility. The expected result (17) 
shows that all these terms are negligible for our purposes. The same holds for Eq. (57). 
We continue and treat the bulk and boundary contributions separately. 



3.3. Bulk contribution 



For the bulk, the quantities IqI^^^^^ {i^ = 1,2) and the function Pi{v) have to be 
calculated. The crucial observation is that Eq. (59), which determines Pi, is well 
known in the study of the spin- 1/2 XXX-Heisenberg chain: Exactly the same function 
determines the T = susceptibility in that model, cf. [10] and references therein. 
Eq. (59) is solved itcrativcly by the Wicncr-Hopf method. The solution reads [16] in 
terms of the Fourier transform V{k) of the function V{v) := P{2uv): 

ft. _i_ fl„ liLB _i_ k = Q 



V{k) = G+{k) 



(61) 



B2' 



with B := B/{2u) and 



iG+(i7r)G^_(0) 
= 16^^^ '^^ 



47r2 



/J,= -%lJ[),^3 = %^(-ln. + l), 



87r3 

The function G+(fc) is given by 



8773 



G+{k) 



(_j^)-ifc/(2,r) 

r(l/2 + ifc/(27r)) 
1 ln(27r) 
2tt 2tt 



-iak 



Combining Eqs. (50), (51) with Eqs. (54) and (58) yields the bulk magnetisation and 
energy in terms of the auxiliary function V{v): 



e — eo = — e 
u 

2-K 

n — no = e 

u 



r(2) 



Q,bulk-'Q,bulk 



[BAO) x(1) 
-'Q.bulk-'Q.bulk 



V{v)e-'"'dv 



V{v)e-''''dv. 



(62) 
(63) 
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Purthermore, from Eq. (38), 



(0) 

Q.bulk 



V{v)dv. 



(64) 



We now successively substitute P{v) from (61) into (64), (62), (63). From the 
substitution of (61) into (64), one obtains 

B= -llnA 

TT So 

G+(i7r) 



So := I, 



(0) 

Q.bulk 



G+(0)- 



(65) 
(66) 



Thus e — eo and n ■ 



no are obtained as functions of s: 

2 



e - eo 



= br 



n — no 



bi 
bs 



u 



:= —I 



•So 

So 

-(0) 



1 + 



bi 



1 + 



(2) 



Ins/so 

bi 
Ins/ So 



+ 62 



In I lns/so| 

In^ s/ So 
In I Ins/sol 



+ 



In s/sq 
b3 



In s/so In s/sq 



(67) 
(68) 



Q,bulk^Q,bulk 
(1) 



U 



2. 



'Q, bulk -'Q, bulk 

:= 2/3i7r/a; 62 := bi/2 
G_(0) 2tt^ 
4 a 

:= G+(i7r)G+(0)/7r 



0s 



2tt 

0J 
a 

'2 

7re 



G^(i7r) 



2/327r' 



, In-i 



a 



The Q-dependence enters only through bQ,CQ. Since we want to calculate the 
susceptibility at constant density, we have to adjust Q such that the density is not 
altered by the finite S-value. This adjustment is done by setting Q = Qq + A, where 
Qo corresponds to B = 00: 

e = eo + (^Qeo) A + 6q/(s/ sq) 
n = no + {dgno) A + CQf{s/so), 
and /(s/so) contains the whole s-dependence of Eqs. (67), (68). From this it follows 



A = 



eo + 



no 
bQ 



f{s/so) 

doeo 



fis/so). 



dqno 

Now the magnetic field h = —dgC and the susceptibility = cf^c- '^'^'^ calculated, 
where x is expressed as a function of h. This calculation is equivalent to considering 
B as a variational parameter and requiring 9b (e — hs) = 0, from which one obtains 
a relation B = B{h). analogous to Eq. (65). This is then substituted into Eq. (64) to 
get s = s{h) and therefrom x = x{h)- This second approach has been chosen in [16]. 
We end up with 

sq ^ ^ 1 In I hih/ho\ 5 

ho 



Xbulk(/l) 

ho 



5Q 



2\nh/ho 
dgeo 



41n^/i//io 161n^/i//io 



dqno 



/so 



(69) 
(70) 
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The constant Sq is given in Eq. (66). Eq. (69) is the key result of this section for 

Xbuik(^)- Note that it is of the form (14) (with E — h there), with specified constants. 
Especially, an alternative expression for Vg has been obtained, cf. Eq. (15): 

^ = 27rvs. (71) 
So 

To complement our discussion of the bulk susceptibility, we are going to express ho in 
terms of the dressed energy functions (35), (36) by combining Eqs. (34), (71). From 
Eqs. (34), (35), (36), it is not difficult to show that at zero magnetic field, 

r'?^e5^^"'=e'(fc)dfc 

2nvs = ^-^ . / . (72) 

/^^e^^^"V(fc)dfc 

An expression for IqI^^^^^ equivalent to Eq. (55) is as follows: 

4:u= £^v(fc)dfc (73) 

rQ 

^{k) =e^ / cosA;'K(i)(sinA;-sinfc')V'(A;')dA;'. (74) 

J-Q 

Thus, by comparing with (35), (36) at /i = 0, 

cQ 



AO) 

-'Q.bulk 



2u [a{v)e^-] = / '=e^(fc) dfc . (75) 

J-Q 

We now insert Eq. (66) into Eq. (71), making use of Eq. (75). Then we obtain 

ho = \— / e^r«'"V,(fc)d/e, (76) 
V eTT J_Q 

which is an expression equivalent to Eq. (70). 



3.4- Boundary contribution 

Let us go back to Eqs. (54), (58). The boundary contribution is calculated by plugging 
these equations into Eqs. (50), (51). The resulting expression is lengthy and we do 
not write it down here. We rather apply the approximation to neglect terms of the 
order 

g-.B/(2„)/52 h 

' In^ s In^ h 

in the ground state energy. These terms would yield a contribution l/(/iln^/i) to 

the susceptibility, cf. Eq. (17). Neglecting these terms means not fixing the scale 

hg in (17) uniquely. However, the algebraic 1/h divergence in (17), dominates such 
terms. Then the boundary contributions read 

[e-eo]B=27r/ a^^v + B)a^^^,^{v + B)dv (77) 
Jo 

POO 

[n-nQ]B = -TT a^,^ii,{v + B)a^^2iki'" + B)dv (78) 
Jo 
1 /"^ 

SB =-- P2{v)iiv, (79) 

4u Jo 
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P2{k) 



where P2 is given by Eq. (60). Note that the only difference with respect to the bulk is 
the different dependence of s on B, cf. Eq. (79). Instead of Sbuik ^ (2xp [—TrB/{2u)] we 
now have Sb ~ 2u/B, which will, according to Eq. (65), cause a logarithmic divergence 
in the magnetisation: 1/B ~ l/ln(/i), yielding the divergence indicated in Eq. (17). 
Terms neglected in (79) arc ^ cxp [—B] ^ ft, + l/(ln/i), and consequently they will 
yield a contribution to the boundary susceptibility which is of the bulk form, Eq. (14). 
These finite terms are negligible compared to the divergent terms given in (17). We 
thus have to insert the expressions for crbuik(^'), '''bulk (^) from the previous section 
into Eqs. (77), (78). The Fourier transform V2{k) of the function 'P2{v) := P{2uv) is 
given in [16], namely 

■ G+(0)(ai/B + a2(lnB)/S2) A; = 
\aiG+{k)/{kB'^), k^O 
1 V2 

Wc now go through the same steps as in the previous subsection: Calculate [e — cq] g , 
[n — no]g, SB as functions of the integration boundary B, and derive therefrom the 
boundary susceptibility as a function of h. This results in 

with the scale ho given in Eq. (70). 
3.5. Explicit expressions in special cases 

We consider three special cases: Half filling for arbitrary coupling, as well as weak and 
strong coupling for arbitrary filling. 

At half fining, g = tt, /g/buik = and 

4'Lk = ^ Tsin^fce^r-Mfc 

Jo 

This expression coincides with the one given by Asakawa et al. [17]. For strong 
coupling, Q = irn and 



ho -\ 

u \ e 



27r3 



^ sin(27rn) 



u — *oo 

Xo ^ — 



^ sin(2'7rn) 



27rn 

1 



=: ho,2 (82) 

(83) 



2?™ 

The latter result has also been obtained in [9]. It is interesting to note that 



ho,2 = Y — "■c, 

where he is the critical field above which the system is fully polarised, Sbuik(/i > he) = 
n/2, [11]. Thus in the strong coupling limit, the logarithmic corrections arc confined 
to fields ft, <C /ic ~ In the special case n = 1 (half filling), Eqs. (82), (83) are 
consistent with known results for the XXX-chain with coupling constant J = 1/u: 
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For this model, x'o^^'^^ = ^/(J''^'^) [10], and the scale /i^^^^^ = A/27r3/e, without 
taking account of the term ~ ln~^ h in (14), [16]. 

The small coupling limit is technically more involved: The integration kernels in 
the integral equations become singular. A numerical evaluation of xoiu) (see section 
3.6 below) confirms the field-theoretical prediction Eq. (A. 15). As far as the scale ho 
is concerned, it is clear that it must diverge for u ^ 0: Exactly at the free-fermion 
point u = 0, the logarithmic corrections in Eq. (14) vanish altogether; corrections to 
Xo at u = arc algebraic with integer powers. The leading contributions ^ /i^, for 
u = are calculated in the Appendix, cf. Eqs. (B.1,B.2,B.5,B.6). 

To describe the divergence of ho\u^o quantitatively, consider first the half- filling 
case where ho ~ e'^°"^*/''. From the general expressions of Iq'^Ii]^, Eq. (55), it is 
clear that this is the case for arbitrary filling. Thus in the small-coupling limit, 
the finite-field susceptibility is obtained from Eq. (14) with |ln/i| <C [In/iol, i-e. 
/lo > /i » l//io ~ exp[-l/u]: 

In/i 

2\nhoJ ' ~ KhFho 

This has to be understood such that the limit u — + is considered at small but finite 
and fixed h. Then the /i-dependent terms are next-leading and can be neglected in a 
first approximation. On the other hand, Xbuik(w — > 0,/i > 0) has been calculated in 
Appendix A, Eq. (A. 14), 

1 2u 

Xbuik(u ^ 0, /i > 0) = + (85) 

2-kvf 1^ vt 



Xbuik(w^ 0,/io > > l//io) = Xo ( 1 + Trm- ) + ^ ( 731- ) -(84) 



with the Fermi velocity Vp = 2sin(7rn/2). Eqs. (84), (85) match provided that 

TT 



ho = const exp 



— Vp 



(86) 



We will confirm numerically this behaviour in section 3.6. 



3.6. Numerical results 

In order to compare the low- field expansion Eqs. (69), (80) to the outcome of the 
non-approximated integral equations (29), (30), we first compute xo, ho numerically. 
To do so, we follow Shiba [9] and rewrite the quantities dqeo, dgno as solutions of 
linear integral equations. Namely, from Eqs. (33), (32), 

rQ 

dqeo = — 4:CosQ p{Q) — 2 / cos kdQp{k)dk 

J-Q 

rQ 

dgno = 2p(Q) + / dQp{k)dk, 
J-Q 

where dQp{k) is obtained from Eqs. (29), (30) at /i = (i.e. B = oo): 

dQp{k) = cos k [^(sin k — sin Q) + K(sin k + sin Q)] p{Q) 
rQ 

+ cos k I K(sin k — sinp) dQp{p) dp. 
J-Q 

These linear integral equations, together with Eqs. (73), (74), are solved numerically. 
The result for xo as a function of u for different fillings is given in Fig. 3, together 
with the small-u expansion (A. 15) and the XXX-limit. We checked numerically 
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Figure 3. The susceptibility Xbulk over u at densities n = 0.2,0.3,0.4,0.5,0.7,1 
(from top to bottom at u = 0). A similar figure has been shown by Shiba [9]. 
Here, we additionally compare with the field-theory result at small u (A. 15) (blue 
dashed lines) and with the XXX-limit (for better comparison the inset shows the 
same figure on a larger scale; the n = 1-line is printed bold here). 



that this way of obtaining xo is equivalent to calculating Vs from Eq. (34) and then 
using Eq. (15). The scale Hq as a function of u at different fiUings is depicted in 
Fig. 4. Besides confirming the small-coupling result Eq. (86), we observe that the 
numerical data are well described by assuming the following form of the constant of 
proportionality in Eq. (86): 

l/n 

7^^^e^ . (87) 

The exponent is exact (cf. Eq. (86)), the constant is conjectured from observing 
good agreement with the numerics, cf. Fig. 4. Having calculated Hq, xo, the next 
step consists in finding x(/i) numerically and comparing with Eqs. (69), (80). The 
calculation of Xbuik(^) is described in [11] (appendix to chapter 6). The idea is 
to rewrite the energy in terms of dressed energy functions (rather than in terms 
of dressed density functions like in Eqs. (29), (30)). The magnetic field enters the 
linear integral equations for the dressed energy functions. Once these equations are 
solved, both the field and the magnetisation are determined. By varying slightly 
the integration boundaries while keeping the density fixed, one performs a numerical 
derivative As/ Ah to obtain x- The results shown in the sequel demonstrate that this 
procedure is highly accurate. 

Fig. 5 shows Xbulk for densities n — 0.2, 1 at different couplings, together with the 
analytical result (69). The boundary contribution X-b(^) can be calculated similarly to 
Xbuik(/i) as sketched above, because the equations are linear in the 1/L contribution. 
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u u 

Figure 4. The scale ho over u for n = 0.2,0.3,0.4,0.5,0.7,1.0 (from bottom to 
top). In a), the large-coupHng result Eq. (82) is verified; in b), the small-coupling 
formula Eq. (87) is visualised (dashed lines). 




Figure 5. The bulk susceptibility for n = 0.2 in a) and for n = 1 in b) at 
u = 1,2,10 (from bottom to top). The dashed blue curves are the analytical 
result (69). The diamonds indicate xOi showing the decrease of the scale with 
increasing u. 
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Results for u = 1, 10 and densities n = 0.2, 1 are shown in Fig. 6. 
4. Friedel oscillations 

Due to the open boundary conditions, translational invariance is broken. This means 
that one-point correlation functions like the density or the magnetisation will no longer 
be just constants but rather become position dependent. In particular, they will 
show characteristic oscillations near the boundaries, the so called Friedel oscillations 
[22]. These oscillations are expected to decay algebraically with distance x from 
the boundary, finally reaching the mean density or magnetisation, respectively, for 
X — > cxD. From a field-theoretical point of view, such a one-point correlation function 
can be obtained from the holomorphic part of the corresponding two-point function 
[34]. We therefore first recall here the asymptotics of two-point functions in the 
Hubbard model according to [35, 11]. After that, we will derive the one-point 
correlation functions for the density n{x) and the magnetisation s^{x). This will 
then allow us to obtain the local susceptibility xi^)- By a conformal mapping we 
will generalise our results to small finite temperatures. As nuclear magnetic resonance 
Knight shift experiments yield direct access to the predictions obtained here 

about its position, temperature as well as density dependence should be valuable for 
experiments on one-dimensional itinerant electron systems. To test our conformal 
field theory results, we present some numerical data based on the density-matrix 
renormalisation group applied to transfer matrices (TMRG). 




Figure 6. The boundary susceptibility XB{h) for u = 1 in a) and for n = 10 in b) 
with n = 0.2, 1 (from top to bottom). The dots are the numerical data, the lines 
the asymptotic form (17) where /Iq^' has been determined by a fit. The insets 
show xg ■ h ■ In^ h, the red horizontal line denotes the asymptotic value 1/4. 
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4-1- Two-point functions 

In the following, we content ourselves with stating the results for the asymptotics of 
pair-correlation functions, without giving any derivations. For any further details, the 
reader is referred to [11, 35]. We also restrict ourselves to the case n ^ 1, i.e., we do 
not consider half-filling. The reason to do so, is that at half-filling the charge sector 
is massive and the correlation functions will be identical to those of the Hcisenborg 
model up to the amplitudes and the spin velocity which do depend on u. The local 
magnetisation and susceptibility for this case, however, have already been discussed 
in [36, 15]. 

The Hubbard model away from half-filling has two critical degrees of freedom 
with diff'erent velocities and the low-energy effective theory outlined in section 2 is 
therefore not Lorcntz-invariant. As the spin and charge excitations arc independent 
from each other we can, however, still apply conformal field theory in this situation 
based on a critical theory which is a product of two Virasoro algebras both with central 
charge c = 1. Due to conformal invariancc, the exponents of the correlation functions 
of primary fields can then be obtained from the finite-size corrections of low-lying 
excitation energies for the finite system. These in turn can be calculated exactly via 
Bethc Ansatz. Then the remaining challenge is to relate the primary fields to the 
original fields of the model. This goal can be achieved by considering the selection 
rules for the form factors involved and by using additional restrictions obtained from 
the Bethe ansatz solution for the finite size spectrum [35]. 

In this situation, the correlation function of two primary fields at zero temperature 
is given by (we include here the dependence on the imaginary time r) 

(<^A±(r,a;)(/.A±(0,0)) = + . .^a- r ^- ^2A+ f ^^iX^ 

with the scaling dimensions 

2A±(AN, D) = f + UDs ± ^E^iMlnisiMl') ' + 2iV± 



2det^ 



2 



2A±(AN, D) = + ^,,D, ± ^ccATV g AiV, > ^ 

\ 2 det 4 / 

Let us explain the symbols used. The entries of the vector AN are integers AA^c, 
ANs, which denote the change of charges and down spins with respect to the ground 
state. The N^^ denote non-negative integers, and D = {Dc, Ds) depends on the parity 
of ANcy. 

Dc = ^ (AA^c + AN^) mod 1 (88) 
= ^AA^c modi. (89) 
Therefore Dc, Dg are either integers are half-odd integers. The matrix ^ has entries 

\ t,sc t,ss J \ ^sc(-D) Z,ss[Q) J 

These entries are obtained from the following system of linear integral equations 
rB 



Zcc{k) = l-\- I Zcs{v)ai{sm.k - v)6.v 

J-B 
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rQ rB 
Zcs{v) = I coskai{v — smk)Zcc{k)dk — / a2{v — v')Zcs{v')dv' . 

J-Q J-B 

Zsc{k) = / ai{smk - v)Zss{v)dv 

J-B 

rQ rB 
■^ss(^) = l+ / coskai{v — s\n.k)Zsc{k)Ak — I a2{v — v')Zss{v')dv' . 

J-Q J-B 

The integration kernels are given by Eq. (31). The integration boundaries B,Q are 
obtained from the linear integral equations for the root densities, Eqs. (29), (30). 

Let us now focus onto (0(r, a;)O(0, 0)) with O = n, .s^, respectively. Since 
the operator O does neither change the particle density nor the magnetisation, we 
have AiVc = = ATV^. Consequently, according to (88), (89), = 0,±1,±2, 
D, = 0,±1,±2,.... Then 

{6iT,x)d{o,o)) - {Of = 

Ai cos{2kF-^x) 

{vcT + ia;)(4o=~?";)^ {vcT — ia;)(C<:c-5<ic)^ (y^j- -\- ia;)(5c<!-5ss)^ (^^7- _ ia;)(io<!-5ss)^ 
^ A2Cos{2kFix) 

(:ori2(/,:j.-: +ki.- ).v 



+ 



{VcT + ix)^cc {VcT — ix)^cc {VsT + \x)^i' {VgT — ix)?" 

A4 COS 2(fci?| + 2kFi)x 

{vcT + ia;)(S<==+€")^ (fcT — ia;)(4<=<^+4''=)^ (zj^t + ixY^'"+i"'>'' {vsT — ia;)(5o»+fs«)^ 

x"^ - v^t"^ , x"^ - U?T^ 
~A^— i^— + Afi 



(91) 



where the amplitudes A^ are different for the density-density and the longitudinal spin- 
spin correlation function. The oscillating terms on the right hand side correspond to 
D = (±l,Tl),(0,±l),(±l,0),(±l,±l) with N ^ {N+,N-,N+,N-) = and the 
non-oscillating to N = (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1) with D = 0. 

4-2. One-point functions in the open system 

From the two-point function {0{zi,Zi)0{z2,Z2)) in Eq. (91) with z = Vc,sT -\- ix we 
can read of the one-point correlation function in the presence of an open boundary by 
considering it as a function of {zi, Z2) only and identifying Z2 = Zi afterwards [34]. 
in(^\\ in(^\\ A cos{2kF^x + (Pi) cos{2kFix + cf>2) 

cos [2(fcf t + fcf -h ^3] . cos [2{kF^ +2kFl)x + (j)4] , , 

^ ' {2x)ii.{2x)iis +^4^2a;)«"+«")=^(2x)(«"+«")=^ ^ ' 

A^ + Ae 
(2a;)2 ' 

with unknown amplitudes Ai and phases (pi where x now denotes the distance from 
the boundary. 

By the usual mapping of the complex plane onto a cylinder we can generalise (92) 
to finite temperatures: 

(d(x)) - (d(x)) -A. cos(2fcfTX + 0i) 
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cos {2kFlX + 4>2) , COS [2(fcFT + ^Fy )x + 



■(^sinh^)«?=(^sinh^)«L ^(^sinh^)«?=(^sinh^)«?. 



+ - 



cos [2(fci?| + 2kF\)x + (}f)4] 
iIi)(«cc+«sc)Vj 

As , ^6 



sinh 2£p)(«cc+«.c)=(^ sinh 2zp)(«„ (^^^ 



(Jk,sinh2i2^)2 (ik, sinh 22:2^)2 ■ 

The magnetic susceptibiHty at zero field is obtained by taking the derivative with 
respect to h ||: Here kp^ = n{n + 2s)/2a, kpi = 7r(n — 2s)/2a and s = xoh for \h\ <C 1 
(without taking account of the logarithmic corrections). Furthermore, the exponents 
and the amplitudes depend on h. However, we neglect this /i-dependence here since 
it gives rise to higher-order contributions in x(a;). In the h = 0-case, 

^cc =: ^, U = l/\/2, Cc. = 0, U = <e/2, 



leading to 



, , „ — ^1 sin(7rnx + + ^2 sin(7rnx + flii2) 
X{x) -Xo = 27rxoa; 



(2a;)fV4(2x)V: 



2 



Aa sm( Snnx + 64) 

+^-xox (;,).w.(2.)V2 • (^4) 

Note that the kp^ + kpi-tcYm in Eq. (92) is /i-independent in lowest order and therefore 
does not contribute to the susceptibility. 

Again we generalise this to finite temperatures: 

-Ai siu(7r;i.i; + Oi) + A2 siii(7r;/,r + (f>2) 



Xix) -Xo= 27rxoa;- 



(^ sinh 2^)«V4(^ sinh 22^)1/2 



„ sin(37rna: + (jjj) 

+ ^^Xox ^.^^ 2^)952/4^^ ^.^^ 2^)1/2 ■ ^^^^ 

Note that we have ignored logarithmic corrections to the algebraic decay of the 
correlation functions throughout this section. Multiplicative logarithmic corrections 
will be present due to the marginal operator in (9). These corrections have been 
discussed for the leading term in (91) in [37]. 



4-3. Numerical results 

To calculate numerically the local magnetisation s^{x) and susceptibility x{x) at finite 
temperatures we use the density-matrix renormalisation group applied to transfer 
matrices (TMRG). The advantage of this method compared to Quantum-Monte- 
Carlo algorithms is that the thermodynamic limit can be performed exactly. This is 
particularly helpful in the present situation where we want to study boundary effects, 
i.e., corrections which are of order 1/L compared to bulk quantities. The method is 
explained in detail in [38, 39, 15]. Here we concentrate on the local magnetisation for 
h^Q and on the local susceptibility for h = 0, both times for generic filling n ^ 1. In 
Fig. 7 TMRG data for {s^{x)) are shown in comparison to the field theory result (93). 
Here the exponents and velocities have been determined exactly by the Bethe ansatz 
solutions Eqs. (90,34) whereas the amplitudes and phases have been used as fitting 

II In the same way the compressibiUty can be obtained by taking derivatives with respect to a chemical 
potential 




Figure 7. TMRG data (black circles) for the local magnetisation (s^(x)) where 
u = 1.0, T = 0.131, s = 0.037 and n = 0.886. In comparison the field theory 
result (93) is shown (blue squares) where the amplitudes and phases have been 
determined by a fit. 



parameters. The agreement is very good. In particular, the exponential decay of the 
correlation function is correctly described by the exponents and velocities obtained by 
Bethe ansatz. The surprisingly rich structure of (s^(x)) is caused by a competition 
between the first three terms in (93) which oscillate with different wave vectors but 
have similar correlation lengths given by — 1.906,^2 — 1.903 and 1^3 = 1.708 (the 
correlation lengths should not be confused with the matrix ^ in (90)). 

In Fig. 8 the local susceptibility for the same set of parameters as in Fig. 7 but 
zero magnetisation is shown. The TMRG data are again compared to the field theory 
result (95) and good agreement is found. Here only the first term of (95) has been 
taken into account because the correlation length belonging to the second term is 
small compared to that of the first one. 

5. Conclusions 

We studied the low-energy thermodynamic and ground-state properties of the one- 
dimensional Hubbard model with open ends. In particular, we concentrated on the 
bulk and boundary parts of the magnetic susceptibility. Based on the low-energy 
effective theory for this model we argued that the functional form of Xhuik, B{h, T = 0) 
and Xbuik,B(^ = 0,T) is universal, i.e., does not depend on filling n or interaction 
strength u. For the bulk susceptibility only the zero temperature and zero field value 
Xo = (27rws)^^ depends on n,u via the spin- wave velocity Vg as does the scale Eq 
appearing in the logarithms. For E <^ Eq with E ^ T,h , however, the scale is not 
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10 20 30 40 

X 

Figure 8. TMRG data (black circles) for the local susceptibility x(a;) where 
u = 1.0, T = 0.131, s = 0.0 and n = 0.886. In comparison the field theory 
result (95) is shown (blue squares) where the amplitudes and phases have been 
determined by a fit. 



important and the bulk susceptibility becomes 

1 / 1 In I In SI 

Xbuik = t; 1 - 



2ttVs \ 2\nE In^ E 
For the boundary part we even find that the result for T = 0, h ^ Hq 

1 / ln|ln/i|\ 

as well as the result for h = Q, T 

1 InllnTl 



12Tlnr V 21nT 

are completely universal. In particular, xb{T — 0,h) and xsiT^h — 0) are identical 
to the results obtained for the Heisenberg model [12, 13]. The same expression for 
Xb{T = 0, h) has also been obtained for the supersymmetric t — J model [18]. The 
universal behaviour of xb at low energies has nothing to do with the special properties 
making the Hubbard model integrable. Instead, the universal behaviour will hold for 
any system whose low-energy effective theory is identical to the one for the Hubbard 
model described in Sec. 2. Therefore even in a generic itinerant electron system, non- 
magnetic impurities or structural defects can give rise to a Curie-like contribution 
to the magnetic susceptibility. This has profound consequences for experiments on 
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such systems, where a Curie term in the susceptibihty is often assumed to be directly 
related to the concentration of magnetic impurities in the sample. In the light of the 
results presented here a more sophisticated analysis is necessary. In particular, the 
temperature dependence of the Curie constant has to be analysed carefully - in the case 
of a boundary considered here the Curie constant will show a logarithmic dependence 
on temperature. In addition, it might be useful to investigate if the Curie contribution 
can be reduced by annealing as one would expect if it is caused by structural defects. 

Based on the Bcthe ansatz solution for the Hubbard model at zero temperature 
we have been able to calculate Xbuik exactly beyond the limit /i <C /lo by determining 
the scale ho for arbitrary filling. In addition, the exact solution has confirmed that 
the bulk and boundary parts show indeed the universal hmctional dependence on 
magnetic field which has been obtained by the low-energy effective theory. 

For the Priedel oscillations in magnetisation and density caused by the open 
boundaries we have derived an asymptotic expansion by making use of conformal 
invariance. We have also calculated the local susceptibility near the boundary 
which is a quantity directly measurable in nuclear magnetic resonance Knight shift 
experiments. We confirmed our results by comparing with numerical data obtained 
by the density-matrix renormalisation group applied to transfer matrices. 
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Appendix A. Small-u expansion 

The Bethe ansatz equations for a finite system are expanded with respect to the 
coupling constant at small couplings, in analogy to the model of an interacting Fermi 
gas [40]. Before turning to the open boundary case, we first perform this expansion for 
periodic boundary conditions. The comparison with open boundary conditions yields 
the surface energy in this approximation. 

Appendix A.l. Periodic boundary conditions 
The energy eigenvalues are given by 

N 

Epbc = -2Y, COS kj (A.l) 
where the kj are obtained through 



= Yi 'r''^'y-'.\ j = i,:;N (A.2) 

^■^ Ai — sm kj — m 

i—i 

^A^-sinfc^^ ff ^'-^ + ;^\ ; = !,.. .,M,. (A.3) 
J-J- A, - sinfc,- - i« -LJ- ^ A, - A„ - 2iM ' ' t V J 

For M = 0, 2M|-many of the are grouped in pairs at 2t:1/L, I — I, ... , M^, and the 
Mj^-many A^'^'s lie at sin(27r//L). The rest of the fcj^^s (namely N — 2M^ many) are 
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not paired, they are at ±2'Kj/L, j = -{N - - l)/2, . . . , -(Af^ + l)/2. We make 
the following ansatz, distinguishing between paired (unpaired) momenta k^^^ (k^^^)- 



k 



(P) _ fc(P.O) 



(A.4) 



(A.5) 



This ansatz is motivated by evaluating numerically the BA equations (A. 2), (A. 3) and 
it is justified a posteriori by observing that the quantities ^j"'^^ , f3j , e; can be obtained 
in a closed form. Eq. (A.4) means that two paired momenta are centred aroTind their 



"centre of mass" k 



(p,o) 



„j ' "J • It can be verified that ei = 6^^^ only for N — 2M| (no 
magnetisation). Expanding the BA equations (A. 2), (A. 3) and comparing coefficients 
of the imaginary and real parts yields 



2^cosfc^°) 



x(0) 



(0) 



Is 



,(0) 



sin fc 



(«,o) 



sin k 



(0) 



cos kj^^ 



t^=-2^cos.rE 



(0) 



sin k: 



(«,o) ■ 



(A.6) 
(A.7) 
(A.8) 



Note that /3j may be positive or negative, depending on the value of cosfc^^"-*. The 
quantity e; in (A.5) may be obtained similarly, however, it will be of no further 
importance for our purposes. Inserting (A.6)- (A.8) into (A.l) results in 



sin cos u , o 

Ej>bc = - 4 ^ ^ + 

sm 5 L 



4.3% 



(A.9) 



where 5 = f - M^. 

Now, with n := N/L, s =: S/L, 



^pbc 



(A.10) 



At 



4 . TTO 
TT T 

(that is, for zero magnetic field), we obtained the charge and spin velocities 



■ sm cos TTS + u (n^ 



45^) 



47r TTO 
^Sm — COSTTS. 



at small u from the low-energy effective Hamiltonian in section 2, 



„ . TTO , 2u 

Vc,s = 2sm — ± — . 

2 TT 



This provides a consistency check on the leading finite-size correction of the ground 
state energy, as obtained from conformal field theory [41] 



^pbc • — 

At finite magnetic field 

h — ^s^pbc 



= e(-) - 



(f c + Vs 



, . TTn . 

-4 sm — sm ws + 2us, 



(A.11) 



(A.12) 



the susceptibility is derived from Eq. (A. 10), 

1 



2u 
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It is important to note that Eq. (A. 13) has been derived at finite s, in the limit of 
vanishing u. In the hmit of smaU fields, s oc h, and thus the smaU- field expansion of 
(A. 13) reads 

1 2u 

Xbuik(w ^0,h>0) = + (A.14) 

with s = s{h). The order of the Umits is important here: The small-coupHng Hmit 
has been taken before the small ft- limit. That is, Eq. (A.14) is valid at small but still 
finite fields, where the field-dependent terms have been neglected (in the main part, 
cf. Eq. (84), a lower bound on the field is given in terms of the scale Hq, namely 
h » l//io)- These singular limits are due to the non-analytic behaviour of Xbuik as a 
function of the magnetic field in the thermodynamic limit, cf. Eq. (17). 

The result for XbuikCi ^ 0, u > 0), that is, with interchanged limits compared to 
Eq. (A.14), is obtained from the low-energy effective Hamiltonian given in section 2: 

Xbulk(/l ^ 0, u > 0) = xo = TT— 

2-KVs 

= + ^ (A-15) 

where the result for Vg given in Eq. (13) has been inserted. The origin of the difference 
between (A.14) and (A. 15) is clarified in section 3.3. 

Appendix A. 2. Open boundary conditions 

For open boundary conditions, the BA equations are given in (24), (25). The remarks 
from Appendix A.l apply similarly to this case, with the modification that all roots 
are on one half of the real axis. Furthermore, the above expansion procedure can be 
repeated with the results 



^1 =-^cosA:(^'°) (A.16) 



^ L + l ^ 



XI A (0) x(0)+oZ]' 



sin k^p 



2(i + l)cosfcf 



(A.17) 



Here the sums run over the symmetrised sets of BA-numbers. We now obtain the 

energy 

. Tr(jV + l) . , ttS 

Eobc = _2!l^ii(i±lL!!l^ + 2 + -^{N^ -4S^ + N- 25). (A.19) 

Expanding (A.19) in powers of l/L yields the boundary contribution to the energy in 
this weak-coupling approximation 

. niT / , . 4 \ ix '^T'' „ 

es = sm — i — 4ssms7r cosstt -|- 2(n — 1) cossttcos — -|- 2 

2 \ TT / 2 

+ u{2s-n){2s + n-l). (A.20) 
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If the boundary susceptibility is derived from this expression as in the previous section, 

one would obtain a constant depending on n and u only. This result cannot be related 
to Eq. (17), showing again the non-commutativity of diverse limits at m ^ 0. 



Appendix B. Free Fermions 



In this section we give Xbuik, Xb for free fermions {u = 0) in the low-energy limit. The 
corresponding quantities are marked by an index 

At T = 0, x''^^\h) is directly obtained from Eqs. (A. 10) (the bulk part) and 
(A. 20), both at u = 0, with the magnetic field given by /i = —dge. In the small-field 
limit, one obtains 

x'^Jlih) = ^ + ^/i^ + 0{h') (B.l) 



1 



167rv|, 



(n - 1) cos ^ 
+ + 



1 7r(n - 1) cos ^ 
2-kvf ' 2vl ' 1&-k\^^ 



+ 0{h^). (B.2) 



Note that XwfkCi) = XB^\h) only for n = l. 

To calculate the susceptibility at finite temperatures, one starts with the free 
energy per lattice site f''^^\ 

L 



j=i L 



In ( 1 -|- exp 



(3 -2 COS' 



TTJ 



1 



H-h/2 



+ {h^ -h) 



,(B.3) 



where the chemical potential /U is to be determined from n ■ 
Euler-MacLaurin formula to Eq. (B.3) yields 

-/?/(//) = (^i + L^i£ In (l + e-/3(-2cosfe-^-/^/2)^ 



-d^x.f^^-^^- Applying the 



1 

"2L 



In (l + e 



0{-2-n-h/2y 



) (l + 6-/5(2-^-^/2)^] ^ ^ (g_4) 



In the T = 0-limit, the results Eq. (B.l), (B.2) are recovered. By performing a saddle- 
point approximation around the two Fermi points in the integral in (B.4), one obtains 
the first T-dependent correction to the zero-field susceptibility. 



x'^^iT) = 
xi'\T) = 



1 



2'kvf 
1 

2-kvf 



+ 



+ 



2-K 



1 - \vl 1 + 



124 



(n — 1) cos 



2 



2vl 



2tt 



rp2 ^ 0(-y4^ 

1 



124 



+ — (n- l)(7 + 3cosn7r) 



cos nTT 7r2(n-l)cos^ 



+ 



34 



(B.5) 



T'^ + 0{T^). (B.6) 



As in the T = 0-case, Xbuik = Xb^' for n = 1 only. 



iff) 
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